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This  study  reports  on  the  development  of  a  FORTRAN  IV  program  to  produce 


double  sampling  acceptance  plans  for  attributes  data.  The  plans  must  satisfy 

cpvt.ix'j  cU*<Lttt»s hi  j  dfka-  fije cfchJt  /f,et 

two  points  on  an  ojer curve,  the-^Ot  point  (p^.  l-^t  )  and  theJJjQt  point 
(P3»  )•  Tw0  models  are  given.  MODEL  I  has  an  additional  constraint  that 

the  maximum  value  of  the  ASN  must  be  less  than  or  equal  to  the  sample  size  for 
a  corresponding  single  sampling  plan.  MODEL  II  relieves  this  constraint.  In 
either  case,  the  resulting  plan  has  a  minimum  ASN  evaluated  at  the  quality 
level  pi  among  all  sampling  plans  satisfying  the  constraints. 


INTRODUCTION 


Frequently  double  sampling  plans  are  employed  In  lieu  of  single  sampling 
plans  for  lot  acceptance  by  attributes  especially  when  lot  or  batch  sizes  are 
large.  A  number  of  systems  of  sampling  plans  are  available  the  most 
recognized  of  which  at  least  in  the  United  States,  are  the  MIL-STD-105D  AOL 


systems  and  the  Dodge-Romig  LTD  and  AOQL  systems. 

In  addition  methods  have  been  reported  for  developing  tailored  double 
sampling  plans  usually  based  on  the  specification  of  two  points  on  an  OC  curve 
(the  likelihood  function).  Two  such  procedures,  largely  taken  from  a  Chemical 
Corps  Engineering  Agency  (1953)  publication,  are  contained  In  Tables  8-2  and 
8-3  of  Duncan  (1974).  These  tables,  based  on  the  Poisson  distribution,  allow 
the  tailoring  of  a  plan  to  a  Producer's  Risk  (a)  of  0.05  at  a  designated 
quality  level  pj  and  a  Consumer's  Risk  of  0.10  at  a  designated  level  p2. 

Plans  may  be  found  for  which  n2,  the  second  sample  size  equals  2nt  and  where 
n2  equals  nj.  In  either  case  the  rejection  number  (cummulative)  on  both 
samples  is  taken  to  be  the  acceptance  number  on  the  second  sample 
(cummulative),  c2,  plus  one. 

In  1969  Guenther,  following  on  some  earlier  work  by  Cameron  (1952), 


developed  what  amounted  to  a  brute-force  algorithm  for  finding  single  sampling 
acceptance  plans  to  satisfy  two  points  on  the  likelihood  function.  The  main 
difference  between  these  two  procedures  is  that  while  Cameron's  procedure 


In  1970  Guenther  extended  his  work  to  double  sampling  plans.  Again  the 


algorithm  is  essentially  brute  force  employing  probability  tables 
extensively.  Unlike  the  Chemical  Corps  tables,  however,  the  fixed 
relationship  between  nj  and  n£  is  not  required  nor  does  the  Poisson  have  to  be 
used.  The  disadvantage  of  his  procedure  is  the  laborious  effort  required  if  a 
plan  is  to  be  developed  by  hand  calculation.  The  algorithms,  however,  are 
sufficiently  simple  to  be  programmed  easily  for  the  computer.  Halley  (1980) 
provides  an  ANSI  Standard  FORTRAN  program  based  on  Guenther's  algorithm  which 
finds  the  minimum  sample  size  single  sampling  plan. 

In  this  study,  a  program  is  developed  for  finding  double  sampling  plans 
and  characteristics  of  the  average  sample  size  functions  of  the  plans  found 
are  compared.  The  basic  algorithm  follows  that  of  Guenther.  However,  an 
objective  function  is  introduced  as  described  in  the  following  paragraphs. 


PROBLEM  FORMULATION 

When  plotted  as  a  function  of  p,  the  likelihood  function  L(p)  of  a 
sampling  plan  provides  the  operating  characteristic  (OC)  curve  for  the  plan. 
Using  the  binomial  distribution  as  an  example,  L(p)  for  a  double  sampling  plan 

is: 


L(p)  = 
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where: 


p  =  incoming  proportion  of  nonconforming  units 
nj  =  first  sample  size 
n 2  =  second  sample  size 
Cj  =  acceptance  number  on  first  sample 
rj  =  rejection  number  on  first  sample 
C2  s  acceptance  number  on  second  sample  (cumulative  for 
"l+n2).  The  rejection  number  on  the  second  sample, 
r2,  is  c2+l. 


In  this  study,  it  is  assumed  that  rj  =  r2  *  c2  +  1.  Thus  the  upper  limit 
of  the  first  summation,  (rj  -  1)  In  the  double  summation  portion  of  L(p)  (the 
probability  of  acceptance  on  the  second  sample),  may  be  replaced  by  c2.  By  so 
doing  a  double  sampling  acceptance  plan  is  fully  described  by  the  plan 
parameters  nlt  n2,  c*,  and  c2. 

Discussion  of  the  algorithm  for  finding  double  sampling  plans  will  break 
L(p)  into  its  two  parts,  the  probability  of  acceptance  on  the  first  sample, 
Pa(nl)’  ar*d  probability  of  acceptance  on  the  second  sample,  Pa(n2).  Thus: 

L(p)  =  Pa(nx)  +  Pa(n2) 
where  for  the  binomial  case: 


Pafn^  =  B(c1  |  nr  p) 


Pa(n2)  =  BB(Cj ,  c^jn^  n?,  p) 
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If  two  design  points  are  designated  on  the  likelihood  function,  ideally  a 
single  (n^,  n0,  Cj,  C2)  set  can  he  found  yielding  an  OC  curve  which  passes 
exactly  through  these  points.  This  would  suggest  setting  the  likelihood 
function  for  each  quality  level  equal  to  its  respective  probabilities  and 
solving  for  a  set  (or  several  sets)  which  exactly  satisfy  the  equations. 
However,  since  nj,  n2*  c^  and  C2  all  must  be  integer,  it  is  doubtful  that  any 
set  can  be  found  which  gives  exact  equality  for  both  equations. 

In  recongnition  of  this  fact  the  Cameron  procedure  for  single  sampling 
plans  selects  a  plan  which  is  as  close  as  possible  to  the  OC  curve  at  the  two 
points.  The  Guenther  procedure  rests  on  the  formulation  of  inequality 
constraints  that  guarantee  risk  levels  at  least  as  good  as  those  stipulated. 

It  is  the  Guenther  procedure  which  is  used  in  this  study. 

The  OC  curve  points  selected  are: 


p^  »  An  Acceptable  Quality  Level  (AQL),  following  the  definition  in 
MIL-STD-105D,  which  should  be  accepted  with  a  probability  of  at 
least  1-a,  a  being  the  Type  I  error  risk. 

P2  «  A  Rejectable  (poor)  Quality  Level  (RQL)  which  should  be  accepted 
with  no  more  than  a  low  probability  B  (Type  II  error  risk). 


These  two  design  parameters  therefore  may  be  expressed  as: 


(plt  l-o)  and  (p2,  6) 


The  resulting  constraint  equations  are: 


L(pt)  >  1-<J 


L(pz)  <  3 


(?) 

(3) 


Actually  an  infinite  number  of  double  sampling  plans  may  be  found  which 
will  satisfy  these  two  inequalities.  Thus  some  measure  of  performance  must  be 
specified  in  order  to  choose  among  them.  The  measure  chosen  in  this  study  was 


to  minimize  the  ASN  when  the  lots  are  at  the  AOL,  pj. 
The  ASN  function  for  a  double  sampling  plan  is: 


ASN  =  n^  +  n2  P(n2) 


where : 

P(n2)  is  the  probability  of  taking  the  second  sample. 
In  binomial  form. 


ASN  =  n.+n,  ^  u 

1  1  dj-c^l  1 
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nx+n2  [8(czln1 ,  p)  -  Btc^n^p)] 


Thus  an  objective  function  was  introduced  as  follows: 

Select  (n^  n?,  Cj.  c2)  to: 

minimize  ASN(pj)  *  nj+n2  [ BCcg  In^pj)  -  B(c1|nr  pj)]  (4) 

Subject  to  equation  (2)  and  (3). 

Two  models  were  formulated  on  the  basis  of  equations  (2),  (3)  and  (4). 
Model  I  introduced  an  additional  constraint,  along  the  lines  of  MIL-STD-1050, 
guaranteeing  that  the  maximum  value  of  the  ASN  of  the  double  sampling  plan, 
ASNMAX,  does  not  exceed  the  sample  size  ns  of  the  minimum  n  and  c  for  a  single 
sampling  plan  satisfying  equations  (2)  and  (3).  This  value  of  n$  is 
designated  ns*. 


In  order  to  study  the  effect  that  the  ASNMAX  constraint  had  on  the  value 

of  the  objective  function,  ASN(pj),  the  constraint  was  removed  in  Model  II  and 

Model  II  was  run  using  the  same  design  parameters.  In  this  way  one  can 

assess,  in  terms  of  average  number  of  items  inspected,  the  penalty  paid  by 

requiring  that  the  ASNMAX  not  exceed  ns*. 

DEVELOPMENT  OF  THE  ALGORITHM 

The  first  step  in  developing  a  double  sampling  plan  is  to  find  the 
minimum  single  sampling  plan  (ns*,  c*)  satisfying  equation  (2)  and  (3) 
where,  for  example: 

Up)  •  l  (>‘l(l-p)"-d  • 

d=0  a 

For  any  fixed  value  of  c,  there  Is  some  maximum  value  of  n  which  will  satisfy 
equation  (2).  If  n  is  increased  beyond  this  value,  designated  nu  for  an  upper 
bound,  equation  (2)  will  be  violated.  Correspondingly,  there  is  a  minimum 
value  of  n,  designated  n^,  which  just  satisfies  equation  (3).  Any  value  of  n 
less  than  n^  will  produce  a  value  of  L(p2)  greater  than  8. 

Beginning  with  c  equals  0,  n^  and  nu  are  found  by  incrementing  n  one  unit 
at  a  time.  If  the  resulting  values  of  ng  is  not  greater  than  or  equal  to 
n^,  c  is  increased  by  one  unit  and  the  procedure  repeated.  Eventually  at 
some  value  of  c,  n,,  >  n,  is  satisfied,  the  current  value  of  c  becomes  c*  and 
n$*  is  set  equal  to  n^.  Thus  the  minimum  single  sampling  plan  satisfying 
equations  (2)  and  (3)  is  found.  ns*  becomes  an  upper  bound  on  nj  for  the 
double  sampling  plan  and  c*  becomes  a  lower  bound  on  c?  where  Cj  must  be  less 


The  logic  behind  these  limits  is  that,  as  Cj  approaches  c*.  nj  will 
approach  ns*.  If  cj  equals  c*.  nj  will  equal  ns*,  03  will  approach  zero  and 
c2  will  equal  Cj.  Thus  any  derived  double  sampling  plan  will  degenerate  to 
the  corresponding  single  sampling  plan. 

The  algorithm  for  deriving  the  double  sampling  plan  employs  equations  (2) 
and  (3)  replacing  the  L(p)  function  with  equation  (1)  or  its  Poisson  or 
hypergeometric  equivalents.  This  study  concentrates  on  the  binomial  of 
equation  (1)  only.  Initially  c2  Is  set  equal  to  c*  and  c^  is  varied  from  0  to 
C2-I.  On  each  iteration  of  Cj,  nj  Is  incremented  from  Cj+1  until  the  equation 

BUjInj,  p2)  <  8 

is  just  satisfied.  This  value  of  n^  becomes  n^,  the  lower  feasible  limit  on 
nl- 

The  algorithm  then  switches  to  the  complete  likelihood  function,  equation 
(1),  searching  to  satisfy 

L(Pj)  >  l-o 

l(p2)  <  &. 

The  computational  procedure  Incorporates  an  Integer  bisection  method  In  order 
increase  efficiency. 

An  upper  bound  on  nj,  nju,  is  obtained  for  the  current  value  of  Cj  from: 
Bfc^n^  p2)  >  1-ct 


with  the  constraint: 


The  algorithm  then  switches  to  the  complete  likelihood  function,  equation 
(1),  searching  to  satisfy: 

n2£:  P2)  ^  BB( c ^ ,c2 | n^ ,n2 *P2)  ^  Q 

n^y :  B(c^ |n^ ,p^)  +  BB(c^ ’C2|n^ ,n2,p^)  ^  1~® 

with  <  rij,  <  n2(J. 

The  lowest  possible  value  for  03  is  ns*  -  nj.  Since  it  is  necessary  to  set  an 
upper  bound  on  n2,  n2u,  in  order  to  use  a  bisection  search  procedure,  an 
arbitrary  value  of  1.5  times  n^  was  used.  This  value  worked  sucessfully  in 
all  cases  examined  herein. 

Once  the  bounds  for  n2  have  been  set,  the  values  of  the  maximum  ASN 
(ASNMAX)  and  of  the  ASN  evaluated  at  pj  (ASN(pj))  are  calcuated  for  each 
feasible  value  of  n2.  These  values  along  with  nj,  n2s  and  n2Jl  form  the  output 
of  the  computer  program. 

ANALYSIS  AND  CONDITIONS  FOR  OPTIMALITY 
The  double  sampling  plan  parameters  which  may  be  varied  are  n^,  n^,  Cj, 
c2,  and  rj,  five  in  all.  Once  minimum  values  for  cj  and  c2  have  been 
established  minimally  satisfying  the  (pj,  1-a)  and  (p2,  B)  points  on  the  OC 
curve,  a  number  of  conbinations  of  nj  and  n2  values  also  will  satisfy  these 
constraints.  Furthermore,  it  is  possible  to  find  feasible  ranges  of  nj  and  n2 
conbinations  for  any  values  of  Cj  and  c2  greater  than  the  minimum  (cj^) 
combination.  Thus  an  infinite  number  of  plans  may  be  found  satisfying  the  two 
OC  curve  points.  They  repesent.  in  effect,  the  solution  of  a  pair  of 
algebraic  equations  in  five  unknowns. 


By  introducing  the  MAXASN  constraint,  the  number  of  feasible  solutions  i 
reduced  as  it  is  if  the  Oodge-Romig  scheme  of  setting  rj  =  r2  *  C2  +1  is 
employed.  It  is  the  introduction  of  the  objective  of  minimizing  the  ASN  at 
the  AQL,  designated  pj,  that  makes  it  possible  to  select  a  single  plan  and 
allow  the  computer  program  to  terminate.  In  order  to  accomplish  this,  the 
performance  of  the  ASN  (p^)  was  analyzed  by  varying  individually  each  of  the 
main  plan  parameters  nj,  n2,  Cj»  and  c2  with  rj  set  equal  to  rg  a  c^  +  1. 

Examples  of  this  procedure  are  presented  in  Appendix  I,  SAMPLE  PROGRAM 
RUNS.  Values  entered  in  each  block  (cell)  are  ASN(pj),  ASNMAX,  nj  and  n2. 
Each  block  represents  the  results  for  a  Cj,  c2  combination  in  the  feasible 
region. 

Some  General  Constraints 

The  value  of  Cj  does  not  exceed  the  c*  of  the  single  sampling  plan  with 
the  same  OC  curve.  If  Cj  equals  c*,  the  double  sampling  plan  will  degenerate 
to  the  single  sampling  plan  with  nj  =  ns*,  n2=0,  and  c2=Cj.  Additionally,  it 
is  obvious  from  the  ASN  formula  that,  if  c^  is  greater  than  c*,  the  ASN(p^) 
will  be  greater  than  ns*  as  will  the  AS WAX.  Hence  the  search  on  c^  may  be 
truncated  at  c*-l. 

Relationship  of  o,  0,  p^,  and  Pg 

For  fixed  a  and  3,  the  acceptance  number  of  a  single  sampling  plan  will 
remain  constant  for  a  constant  discrimination  ratio,  0,  defined  as  the  ratio 
p2/pl.  As  pi  increases,  the  sample  size  of  a  single  sampling  plan  decreases 
considerably.  This  can  be  explained  by  the  fact  that  it  is  the  absolute 
difference  (p2-pl)  that  influences  n,  not  the  discrimination  ratio.  For 
example,  for  pl=0.02  and  p2=0.10  the  single  sample  size  is  approximately  38 


with  an  acceptance  number  of  3.  For  pl=0.001  and  p2*0.005,  the  same 
discrimination  ratio,  the  single  sample  size  is  approximately  1,350  with  an 
acceptance  number  of  3.  This  confirms  that  if  (p2-pl)  is  small,  n  will  be 
large  and  if  (p2-pl)  is  large  n  will  be  small.  Furthermore  as  the 
discrimination  ratio  decreases  to  1.5,  the  acceptance  number  increases  rather 
dramatically  to  52.  This  result  is  verified  by  use  of  the  Poisson 
approximation  to  the  binomial  and  the  methodology  used  to  solve  for  single 
sampling  plans  therewith.  See  Cameron  (1952). 

Larger  values  of  cl  result  in  a  smaller  range  between  nls  and  nil.  As  cl 
increases  for  a  specific  value  of  c2,  the  lower  bound  on  nl  Increases.  This 
in  turn  reduces  the  number  of  double  sampling  plans  computed  in  each  cell 
because  the  upper  bound  on  nl,  nlu,  is  dependent  on  c2,  not  on  cl. 

Effect  of  ASNMAX  Constraint,  MODEL  I 

When  the  constraint  ASNMAX  <ns*  is  Imposed,  the  values  of  ASN(pl)  In  the 
region  where  c2  exceeds  c*  become  Infeasible.  An  advantage  of  this  property 
is  that  the  search  routine  to  locate  the  global  minimum  does  not  need  to 
search  the  region  where  c2  exceeds  c*.  However,  MODEL  II,  In  which  the  ASNMAX 
constraint  is  not  imposed,  requires  the  evaluation  of  columns  for  c2  >  c*.  It 
is  practically  important  to  study  the  difference  in  minimum  ASN(pl)  in  each 
case  until  a  global  minimum  has  been  found.  Thus,  MODEL  II  searches  for 
column  minimums  and  selects  the  global  minimum  from  that  group.  MODEL  I  needs 
only  to  look  at  the  values  for  c2=c*. 

Behavior  of  ASN(p^)  as  a  Function  of  n^ 

Plotting  of  the  values  of  ASN(p|)  as  a  function  of  nj  showed  that  it  is  a 
quasi-convex  function  of  n^  ( Integerization  of  nj  and  n?  may  explain  why  the 


results  were  not  purely  convex)  Thus  a  search  for  the  minimum  ASN(pj)  needs 
only  to  continue  one  step  beyond  the  point  at  which  the  minimum  exists. 


Behavior  of  ASN  (p^)  as  a  Function  of  ct  . 

In  general,  the  ASN(pj)  proved  to  be  a  quasi-convex  function  of  Cj. 

For  a  constant  discrimination  ratio,  (D  =  p2/pl),  MODEL  I  results  showed 
that  the  minimum  ASN(Pl)  occurs  at  the  same  value  of  cl  irrespective  of  the 
value  of  pi.  Such  is  not  the  case  with  MODEL  II. 

The  values  of  Cj  where  the  minimum  ASN(pl)  occurs  is  influenced  by  the 
ratio  p2/pl.  As  the  ratio  decreases,  the  value  of  cl  Increases.  However  this 
relationship  also  is  affected  by  the  magnitude  of  (p2-pl). 


Table  I. 


DESCRIM I NATION 
RATIO  (P2/P1) 

VALUE  of  cl 

AT  min.ASN(pl) 
of  column.(c2*c*) 

25 

0 

10 

0  or  1 

5 

0  or  1 

4 

1 

3 

2 

2.5 

2  or  3 

Table  I  indicates  that  as  the  ratio,  D,  decrease,  the  minimum  ASN(pl) 
tends  to  increase.  That  is,  the  corresponding  value  of  cl  becomes  larger.  As 
the  ratio  Increases  curves  of  ASN(pj)  shift  and  become  truncated  constantly 
increasing  from  the  first  feasible  solution  rather  than  moving  downward  to  a 
minimum  before  sweeping  upward.  In  MODEL  II,  as  c2  Increases,  the  value  cl 
for  the  minimum  ASN(pl)  of  each  column  may  vary. 


Behavior  of  ASN(pj)  as  a  Function  of  c2 

Plotting  of  the  ASN(p^)  values  as  a  function  of  c2  yielded  quasi-convex 
results  as  well.  However,  there  were  substantially  different  results  under 
the  two  models.  When  MODEL  I  (constrained  ASNMAX)  was  employed,  the  minimum 
ASN(pl)  value  occured  always  for  values  of  c2  equal  to  c*.  Under  MODEL  II, 
this  was  the  case  occassional ly  but  not  always. 

Impact  of  the  ASNMAX  Constraint 

The  main  objective  of  the  analysis  of  MODEL  II  was  to  evaluate  the  effect 
of  the  constraint  ASNMAX<ns*  on  the  objective  function. 

Table  II  shows  the  minimum  ASN(pl)  obtained  from  both  models  together 
with  their  ASNMAX  values  for  some  representative  cases. 


TABLE  II 


VALUE 

OF 

Pi 

VALUE 

OF 
—  P2 

rat  nr 

OF 

P2/Pl 

“MooEU'r 
ASN(Pl)  ASN 

MAX 

“MODEUTI 
ASN(Pl)  I  ASN 

1  max 

- ? - 

OF 

REDUCTION 

0.001 

0.04 

A 

179.8 

185. 

149.1 

232. 

17.1 

0.02 

0.8 

4 

75.4 

86.5 

73.2 

115. 

2.23 

0.005 

0.125 

25 

19.9 

27.1 

19.9 

27.1 

0 

0.005 

0.05 

67.1 

98.9 

67.1 

98.9 

0 

0.04 

0.2 

5 

21.6 

31.5 

21.6 

31.5 

0 

As  indicated,  a  reduction  in  ASN(pl)  generally  is  obtained  only  when  the  D 
ratio  is  very  low  and/or  the  difference  between  p2  and  pi  is  small.  For  a  D 
ratio  of  4  and  (p2-pl)  equals  0.03,  a  reduction  in  minimum  ASN(pl)  of  17.1%  is 
obtained  by  dropping  the  constraint.  However,  when  the  difference  between  p2 
and  pi  is  increased  to  0.6  with  the  same  D,  a  reduction  of  only  2.9%  in 


ASN(pl)  is  seen,  but  seen  at  a  cost  of  a  substantial  increase  in  ASNMAX.  In 
the  rest  of  the  cases  illustrated,  no  reduction  in  minimum  ASN(pl)  is  obtained 
by  eliminating  the  ASNMAX  constraint.  For  these  cases  the  D  ratio  was  5  or 
greater  and  (p2-pl)  was  0.045  or  greater.  These  results  indicate  that  both 
the  0  ratio  and  the  difference  (p2-pl)  affect  ASNMAX  but  only  when  both  are 
small . 

Additionally,  whenever  a  reduction  in  minimum  ASN(pl)  is  obtained  by 
dropping  the  constraint,  the  Increase  in  corresponding  ASlfiAX  value  may  be 
substantial.  However  the  Increase  In  ASNMAX  may  become  smaller  as  the 
difference  between  p2  and  pi  increases.  In  other  words,  as  the  difference 
between  p2  and  pi  decreases,  the  price  to  be  paid  for  the  protection  against 
high  ASNMAX' s  will  start  to  increase. 

SUMMARY  OF  THE  ALGORITHMS 

MODEL  I 

STEP  1.  Compute  smallest  single  sampling  plan,  (ns*,  c*). 

STEP  2.  Set  c2  =  c*. 

STEP  3.  Incrementing  on  cl(0,  1,  2,...,  c*-l): 

3a.  Compute  feasible  bounds  on  nl;  i.e.,  nil  and  nlu. 

3b.  Incrementing  on  nl(nll,  nll+1,...,  nlu)  compute  bounds  on  n2; 
i.e.,  n21  and  n2u. 

3c.  Incrementing  on  n2(n21,  n21+l,...,n2u)  compute  ASNMAX  and  ASN(pl). 

Condition  for  Optimality 

Feasible  values  are  those  for  which  the  likelihood  constraints,  (pl,l-a) 


and  (p2,3),  and  the  ASNMAX  constraints  are  satisfied.  At  each  calculation  in 
the  feasible  region,  ASN(pl)  is  calculated  and  the  optimal  double  sampling 


plan  is  the  one  with  min  ASN(pl).  Because  of  convexity  of  ASN(pl),  the 
algorithm  shifts  from  cell  to  cell  (value  of  Cj)  whenever  the  current 
calculation  (ASN(pl)'s)  exceeds  that  for  the  previous  calculation. 

MODEL  II 

STEP  1.  Compute  the  smallest  single  sampling  plan  (ns*,  c*). 

STEP  2  Incrementing  on  c2(c*,  c*+l,  c*+2,...): 

STEP  3  Then  incrementing  on  cl(0,  1,  2,,,,c*-l): 

3a.  Compute  feasible  bounds  on  nl;  i.e;  nil  and  nlu. 

3b.  Incrementing  on  nl(nll,  nll+1 .  nlu)  compute  bounds  on  n2; 

i.e;  n21  and  n2u. 

3c.  Incrementing  on  n2(n21,  n21+l,  n21+2, . . .  ,n2u)  compute  ASN(pl). 

Condition  for  Optimality 

Feasible  values  are  those  for  which  the  likelihood  constraints  are 
satisfied.  At  each  calculation  in  the  feasible  region,  ASN(pl)  Is  calculated 
and  the  optimal  double  sampling  plan  is  the  one  with  min  ASN(pl).  Because  of 
convexity  of  ASN(pl)  the  algorithm  shifts  from  cell  to  cell  (value  of  Cj) 
until  the  current  calculation  (ASN(pl)'s)  exceeds  that  for  the  previous  cell 
calculation.  Slmllarlly  the  algorithm  shifts  column  to  column  (value  of  C2) 
until  the  minimum  ASN(pl)  of  the  current  column  exceeds  that  for  the  previous 
col umn . 


COMPUTER  CODE 

The  program  originally  was  written  in  FORTRAN  IV  to  run  on  a  PDP  11-34 
computer.  Later  It  was  modified  to  allow  r^  to  be  entered  externally  (rather 
than  set  at  ^  *  C£+l)  and  to  run  on  a  VAX  11-750  computer.  The  complete  code 
is  Included  in  APPENDIX  II. 


As  stated  previously,  the  single  sampling  plan  is  computed  using  a  brute 
force  method;  i.e:  the  search  starts  with  an  acceptance  number  of  zero  and  the 
sample  size  is  Incremented  by  one  at  each  iteration  until  L(pl)  and  L(p2) 
satisfy  their  respective  inequalities.  If  the  solved  value  of  nl  exceeds  nu, 
no  feasible  solution  exists  for  that  value  of  c,  c  is  incremented  by  one,  and 
the  search  process  for  nl  and  nu  starts  anew.  Depending  on  the  input 
parameters  (a,  $,  pi  &  p2),  the  single  sample  size  may  become  very  large  thus 
requiring  a  large  number  of  iterations  to  reach  the  first  feasible  solution. 

To  reduce  unnecessary  computations,  the  user  may  input  a  "seed"  number  as  a 
starting  value  of  the  single  sample  size.  The  closer  the  seed  is  to  the  true 
solution,  the  lesser  the  number  of  Iterations  required.  However,  the  user 
must  be  very  careful  in  entering  a  seed  value.  If  a  higher  value  of  the  seed 
than  the  true  solution  minimum  ns  is  entered,  the  algorithm  will  converge  to  a 
single  sampling  plan  with  an  acceptance  number  higher  than  that  of  minimum 
single  sampling  plan  (the  desired  solution). 
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APPENOIX  I 
SAMPLE  PROGRAM  RUNS 


DEPT.  OF  ISE 
UNIVERSITY  OF  FLORIDA 


*****DOUBLF  SAMPLING  SYSTEM***** 


ALPHA  =0.  0500  BETA  =0.  1000 

PO  =0.0100  PI  =0.0400 

REJECTION  NO.  OF  FIRST  SAMPLE  (Rl)  = 


S -  ■ 

ACCEPTANCE  NO.  <C)  »  4 
LOWER  BOUND  ON  N  (NS)  =  198 
UPPER  BOUND  ON  N  (NL)  =  198 

DOUBLE  SAMPLING  PLANS 


FOR 

r> 

ll 

o 

C2= 

4 

136 

62 

62 

186.  8584 

137 

61 

61 

187.  0367 

DOUBLE  SAMPLING 

PLANS 

FOR 

Cl«  1 

C2= 

4 

164 

34 

34 

185.  1312 

165 

33 

33 

185.  5088 

DOUBLE  SAMPLING 

PLANS 

FOR 

Cl=  2 

C2= 

4 

183 

15 

IS 

189.  1240 

184 

14 

14 

189.  7155 

DOUBLE  SAMPLING 

PLANS 

FOR 

Cl=  3 

C2= 

4 

188 

11 

11 

190.  1723 

194 

4 

4 

194.  7896 

C2-M  1> 


181.  4362 
181.  8368 


179.  7643 

180.  3866 


186.  5909 

187.  3789 


188.  8725 
194.  3390 


i  Pitmen  r.Tj  i> 


ACCEPTANCE  NO.  (C)  »  5 
LOWER  BOUND  ON  N  (NS) 
UPPER  BOUND  ON  N  (NL) 

DOUBLE  SAMPLING  PLANS 


FOR 

C I  =  0 

C2=  5 

63 

213 

213 

251.  3808 

162.  9084 

64 

207 

212 

247.  0489 

162.  1922 

65 

203 

210 

244.  4871 

162.  3609 

DOUBLE  SAMPLING  PLANS 


FOR 

C 1  =  1 

C2=  5 

101 

180 

183 

232.  8200 

149.  1270 

102 

174 

181 

229.  4124 

149.  1608 

DOUBLE  SAMPLING  PLANS 

FOR 

Cl-  2 

C2-  5 

134 

161 

162 

223.  2362 

157  9882 

135 

151 

160 

218.  6858 

157.  8415 

136 

143 

158 

215.  2449 

157  9577 

DOUBLE  SAMPLING  PLANS 

FOR 

Cl=  3 

C2=  5 

166 

149 

155 

220.  4352 

177.  8430 

167 

127 

151 

213.  3944 

177.  2542 

168 

114 

148 

209.  6420 

177  3490 

DOUBLE  SAMPLING  PLANS 
FOR  C 1 =  4  C2—  5 

198  134  262  221.  8191 

199  87  262  214.  4609 

200  72  262  212. 7945 

GLOBAL  MINIMUM  ASN(P0>«  149.  16 

CORRESPONDING  N1  =  102 

CORRESPONDING  N2S  -  174 

CORRESPONDING  Cl  *  1 

CORRESPONDING  C2  =5 


202.  6424 
202.  0609 
202.  5721 


DEPT.  OF  ISE 
UNIVERSITY  OF  FLORIDA 


*****DOUBLE  SAMPLING  SYSTEM***** 

ALPHA  =0.  0500  BETA  =0.  1000 

PO  =0.  0200  PI  =0.  0800 

REJECTION  NO.  OF  FIRST  SAMPLE  (Rl)  =  C2+(  1) 

S - - - 

ACCEPTANCE  NO.  (C>  =  4 
LOWER  BOUND  ON  N  (NS)  -  98 
UPPER  BOUND  ON  N  <NL)  “  99 

DOUBLE  SAMPLING  PLANS 
FOR  Cl*  0  C2*  4 

46  54  54  90.  6816  78.  5600 

47  53  53  90.  8418  79.  3637 

DOUBLE  SAMPLING  PLANS 
FOR  Cl*  1  C2*  4 

62  39  39  86.  5216  75.  4389 

63  38  38  86.  8856  76.  3481 

DOUBLE  SAMPLING  PLANS 
FOR  Cl*  2  C2=  4 

75  27  27  86.  1676  79.  6567 

77  24  24  86.  9210  81.  3386 

DOUBLE  SAMPLING  PLANS 
FOR  Cl*  3  C2=  4 

87  16  16  90.2002  88.0654 

88  14  15  90.7994  88.9571 


— -  <  mm..  w  mm  MMjh'C'yp 


.-jk  rJ^j.^jkTVT'jnarjr-j.:  _ . 


S - 

ACCEPTANCE  NO.  <C)  «  5 
LOWER  BOUND  ON  N  <NS>  =  114 
UPPER  BOUND  ON  N  (NL)  =  131 

DOUBLE  SAMPLING  PLANS 


FOR  C 1 =  O  C2=  5 


31 

105 

108 

124.  7199 

79.  8665 

32 

100 

106 

121.  2064 

79.  6080 

33 

96 

105 

1 18.  5924 

79.  7088 

DOUBLE  SAMPLING 

PLANS 

FOR 

Cl=  1 

C2= 

5 

50 

88 

93 

115.  1913 

73.  2101 

51 

82 

91 

111.  7191 

73  2325 

DOUBLE  SAMPLING 

PLANS 

FOR 

Cl=  2 

C2= 

5 

66 

81 

84 

111.  4853 

77.  6491 

67 

72 

81 

107.  4158 

77.  6818 

DOUBLE  SAMPLING 

PLANS 

FOR 

Cl=  3 

C2= 

5 

82 

75 

83 

109.  7731 

87  7234 

83 

57 

79 

104.  1007 

87.  4925 

84 

49 

75 

102.  1335 

87.  9863 

DOUBLE  SAMPLING 

PLANS 

FOR 

Cl=  4 

C2= 

5 

98 

63 

131 

109.  3472 

100.  0915 

99 

35 

131 

105.  3023 

100.  1993 

GLOBAL  MINIMUM 

ASN(PO) 

=  73.  23 

CORRESPONDING  N1  =51 

CORRESPONDING  N2S  =  82 

CORRESPONDING  Cl  =1 

CORRESPONDING  C2  =5 


20 


,v, 


DEPT.  OF  ISE 
UNIVERSITY  OF  FLORIDA 


*****DOUBLE  SAMPLING  SYSTEM***** 


ALPHA  =0.  0500 
PO  *0.  0200 


BETA  =0.  1000 
P 1  =0.  0600 


REJECTION  NO.  OF  FIRST  SAMPLE  (Rl)  =  C2+< 


ACCEPTANCE  NO.  (C)  =  7 
LOWER  BOUND  ON  N  (NS) 
UPPER  BOUND  ON  N  (NL) 


DOUBLE  SAMPLING  PLANS 


FOR  C 1=  O  C2=  7 


FOR  Cl*  1  C2=  7 


DOUBLE  SAMPLING  PLANS 


Cl*  2  C2=  7 


DOUBLE  SAMPLING  PLANS 


FOR  Cl> 


DOUBLE  SAMPLING  PLANS 


FOR  Cl*  4  C2=  7 


DOUBLE  SAMPLING  PLANS 


FOR  Cl*  5  C2=  7 


DOUBLE  SAMPLING  PLANS 


FOR  Cl*  6  C2=  7 


195.  7031 
194.  7807 


DOUBLE  SAMPLING  PLANS 


187.  7091 
186.  9561 


180.  3470 
179.  8239 


174.  2581 
174.  0120 


172.  2300 
172.  2892 


174.  6661 
174.  4244 


181.  5391 
181.  4739 


152  9495 

153  6086 


139.  1258 
139.  9612 


135.  4178 

136.  3186 


140.  0145 
140.  9156 


148.  9510 

149.  8648 


161.  6958 

162.  4568 


I  76.  6320 
177.  3945 


ACCEPTANCE  NO  (C)  =  B 
LUWtK  UUUND  UN  N  t  NS )  =  in  5 
UPPER  BOUND  ON  N  (NL)  =  236 

DOUBLE  SAMPLING  PLANS 
FOR  Cl=  0  C2=  8 

44  197  199  234.  9933  160  0134 

45  193  198  232.  0821  160  2446 

DOUBLE  SAMPLING  PLANS 
FOR  C 1 =  1  C2=  8 

70  173  176  227.3926  140.8520 

71  169  174  224.7274  141  38^5 

DOUBLE  SAMPLING  PLANS 
FOR  Cl=  2  C2=  8 

92  156  157  220.2482  135.6124 

93  151  156  217.1164  136.0340 

DOUBLE  SAMPLING  PLANS 
FOR  Cl*  3  C2=  8 

113  140  141  212.4149  139.6817 

114  134  139  209.1393  140.0779 

DOUBLE  SAMPLING  PLANS 
FOR  Cl*  4  C2*  8 

134  120  125  203.  5589  149.  6677 

135  113  123  200.4914  1500802 

DOUBLE  SAMPLING  PLANS 
FOR  Cl*  5  C2=  8 

154  106  116  200.3735  163.1247 

155  95  113  196.5552  163.3638 

DOUBLE  SAMPLING  PLANS 
FOR  Cl*  6  C2=  8 

174  91  125  200.3812  178.8247 

175  75  11B  196.7400  179.0665 

DOUBLE  SAMPLING  PLANS 
FOR  Cl*  7  C2*  8 

194  76  236  204.  8341  195  9521 

195  49  236  201.9844  1962862 


zt 


M 


a 


ACCEPTANCE  NO.  <C)  «■  9 
LOWER  BOUND  ON  N  (NS) 
UPPER  BOUND  ON  N  (NL> 

DOUBLE  SAMPLING  PLANS 


FOR 

O 

II 

H 

u 

C2=  9 

41 

236 

242 

272.  4937 

173.  9185 

42 

229 

240 

266.  3910 

1/2  9/66 

43 

224 

239 

262.  6524 

173.  0344 

DOUBLE  SAMPLING  PLANS 


FOR 

C 1  =  1 

C2=  9 

66 

221 

222 

273. 3678 

150.  2853 

67 

212 

220 

265.  8895 

149.  3582 

68 

205 

218 

260.  2918 

149.  0869 

69 

200 

217 

256  5717 

149  5144 

DOUBLE  SAMPLING  PLANS 


FOR 

C 1  =  2 

C2=  9 

89 

204 

204 

266.  7397 

142.  7357 

90 

193 

202 

258.  1279 

141.  8807 

91 

185 

200 

252.  1336 

141.  7317 

92 

179 

198 

247.  8829 

142.  0573 

DOUBLE  SAMPLING  PLANS 


FOR 

Cl«  3 

C2=  9 

111 

186 

188 

256.  2776 

145  0258 

112 

173 

186 

247.  1019 

144.  3393 

113 

164 

184 

241.  0539 

144.  3170 

114 

157 

182 

236.  5696 

144.  6165 

DOUBLE  SAMPLING  PLANS 

FOR 

C 1=  4 

C2=  9 

132 

175 

176 

249.  4611 

154.  0411 

133 

156 

174 

237.  6933 

153.  0990 

134 

145 

171 

231.  2966 

153.  1052 

DOUBLE  SAMPLING  PLANS 

FOR 

Cl=  5 

C2=  9 

153 

154 

171 

237.  1488 

166.  3868 

154 

132 

167 

226.  1181 

165.  7449 

155 

121 

164 

221.  0988 

166.  0171 

',rrvi 


DOUBLE  SAMPLING  PLANS 
FOR  Cl=  6  C2=  9 

174  120  189  223.  4432  181.  0777 

175  103  181  217.  4335  181.  2191 

DOUBLE  SAMPLING  PLANS 


FOR 

Cl*  7 

C2=  9 

194 

112 

273 

224.  6025 

198.  0900 

195 

81 

273 

217.  1298 

198.  0276 

196 

70 

273 

215.  1222 

198.  6775 

DOUBLE  SAMPLING  PLANS 
FOR  Cl*  8  C2*  9 

215  53  273  222.  1338  215.  9652 

216  41  273  221.5181  216.7635 

GLOBAL  MINIMUM  ASN(PO)*  135.42 

CORRESPONDING  N1  -  101 

CORRESPONDING  N2S  *  105 

CORRESPONDING  Cl  =2 

CORRESPONDING  C2  =7 


24 


DEPT.  OF  ISE 
UNIVERSITY  OF  FLORIDA 

»****DOUBLE  SAMPLING  SYSTEM***** 

ALPHA  =0.  0500  BETA  =0.  1000 

PO  =0.0150  PI  =0.0450 

REJECTION  NO.  OF  FIRST  SAMPLE  (Rl)  =  C2+<  1) 

S - - - 

ACCEPTANCE  NO.  (C)  =  7 
LOWER  BOUND  ON  N  <NS>  =  260 
UPPER  BOUND  ON  N  (NL)  =  266 

DOUBLE  SAMPLING  PLANS 
FOR  Cl=  O  C2=  7 

75  195  195  260.  0430  207.  2268 

76  193  194  259.  1314  207.  8022 

DOUBLE  SAMPLING  PLANS 
FOR  Cl=  1  C2=  7 

108  164  164  250.  0084  187.  1743 

109  162  162  249.  2654  187.  9875 

DOUBLE  SAMPLING  PLANS 
FOR  Cl=  2  C2=  7 

138  135  135  239.  5929  184.  0606 

139  133  134  239.  0794  184.  9139 

DOUBLE  SAMPLING  PLANS 
FOR  Cl-  3  C2=  7 

164  111  111  232.4163  189.4726 

165  109  109  232.1787  190.3486 

DOUBLE  SAMPLING  PLANS 
FOR  Cl-  4  C2—  7 

188  90  90  229.  9069  201.  1494 

189  87  88  229.5068  201.8962 

DOUBLE  SAMPLING  PLANS 
FOR  Cl*  5  C2=  7 

211  72  72  233.  1982  217.  1050 

212  68  70  232.9633  217.8536 

DOUBLE  SAMPLING  PLANS 
FOR  Cl=  6  C2=  7 

234  59  64  242.  9257  236.  2712 

235  51  60  242.7148  236  9932 


1, 


iWWiWJ 


IWCW 


‘lV»V  I1, 


.1 


s - 

ACCEPTANCE  NO.  (C>  -  9 
LOWER  BOUND  ON  N  (NS)  =  314 
UPPER  BOUND  ON  N  <NL>  -  363 

DOUBLE  SAMPLING  PLANS 


FOR 

Cl-  0 

C2-  9 

55 

316 

322 

56 

309 

320 

57 

303 

319 

58 

298 

317 

DOUBLE  SAMPLING  PLANS 


FOR 

Cl=  1 

C2=  9 

89 

291 

294 

90 

283 

292 

91 

276 

290 

92 

270 

289 

DOUBLE  SAMPLING  PLANS 


FOR  Cl-  2  C2=  9 


120  264  269 

121  255  267 

122  247  265 

123  241  263 

DOUBLE  SAMPLING  PLANS 


FOR 

Cl-  3 

C2—  9 

149 

243 

248 

150 

231 

246 

151 

222 

244 

152 

215 

242 

DOUBLE  SAMPLING  PLANS 

FOR 

Cl-  4 

C2—  9 

177 

226 

232 

178 

209 

229 

179 

197 

227 

180 

189 

224 

DOUBLE  SAMPLING  PLANS 

FOR 

Cl-  5 

C2=  9 

205 

198 

223 

206 

178 

220 

207 

166 

216 

DOUBLE  SAMPLING  PLANS 

FOR 

Cl-  6 

C2-  9 

232 

200 

248 

233 

156 

241 

234 

139 

233 

364  4188  233  3808 
358.  5363  232.  4478 
353.6355  231.9713 
349.7147  231  9728 


361.2587  201.4461 

354.  7504  200.  8535 

349.  1775  200.  5675 

344.  5428  200.  6042 


349.  1373  191  0035 

342.  3050  190.  6154 

336.  3438  190.  4330 

332.1181  190.7493 


337.  9465  194.  2755 

329.  6002  193.  7339 

323.  5883  193.  7005 

319.  1306  194.  0067 


327.  9336  206.  0636 

317.  5687  205.  3329 

310.  5440  205.  1960 

306.  1929  205.  5502 


312.  6136  222.  6100 

302.  7364  222.  1069 

297.  2081  222.  2803 


313.  9585  243.  8767 
296.  9233  242  4273 
290  9528  242  5471 


DOUBLE  SAMPLING  PLANS 


P 


SR 

y. 

y. 

jft 

5 


FOR  Cl-  7  C2-  9 


260 

261 


126 

106 


363 

363 


DOUBLE  SAMPLING  PLANS 


FOR  Cl-  8  C2-  9 


294.  2362 
289.  8002 


264.  7458 

265.  0618 


287  82  363  297.  9774 

288  63  363  296.  4333 

GLOBAL  MINIMUM  ASN(PO)-  184.06 

CORRESPONDING  N1  =  138 

CORRESPONDING  N2S  -  135 

CORRESPONDING  Cl  -  2 

CORRESPONDING  C2  -  7 


288.  5138 

289.  1825 
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APPENDIX  II 
FORTRAN  IV  PROGRAM 


0001 

0002 

ooo:j 

0004 

0005 

0006 

0007 

0008 

0009 

0010 

OOll 

0012 

0013 

0014 

0015 

0016 

0017 

0018 

0019 

0020 

0021 

0022 

0023 

0024 

0025 

0026 

0027 

0028 

0029 

0030 

0031 

0032 

0033 

0034 

0035 

0036 

0037 

0038 

0039 

0040 

0041 

0042 

0043 

0044 

0045 

0046 

0047 

0048 

0049 

0050 

0051 

0052 

0053 

0054 

0055 

0056 

0057 


QUALITY  CONTROL  DOUBLE  SAMPLING  PROGRAM  TO  ANALYSE 
DOUBLE  SAMPLING  PLANS.  ASN(PO)  AND  ASNMAX. 

BINOMIAL  AND  POISON  PROBALITY  DISTRIBUTIONS  USED. 

PROGRAMED  BY  R.  WAREN  RANGARAJAN 
INDUSTRIAL  AND  SYSTEMS  ENGINEERING  DEPARTMENT 
UNIVERSITY  OP  FLORIDA 
GAINESVILLE.  FLORIDA  32611 

DOUBLE  PRECISION  SUMLOG 

INTEGER  C.  Cl.  C2,  C1MIN.  C2MIN.  Rl.  Rll 

BYTE  STING<8> 

C0MM0N/BLK1/N2S, N2L 
C0MM0N/BLK2/PS. PL 
C0MM0N/BLK3/N 1 
C0MMQN/BLK4/ ALPHA,  BETA 
C0MM0N/BLK5/P0,  PI 
C0MM0N/BLK6/C 1 .  C2 
C0MM0N/BLK7/SUML0G ( 2500 ) 

COMMON /BLK8/N 

C0MM0N/BLK9/C2MAX .  C1MAX( 15) 

COMMON/BLK 10/NS,  NL 
C0MM0N/BLK1 1 / ASN.  ASNMAX 

WRITE (5. *)  '  NAME  OF  OUTPUT  FILE?' 

READ (5, 1)  STING 
1  FORMAT ( 10A1 ) 

CALL  ASSIGN  (1. STING) 

BEGINNING  INITIALIZATION 

N-0 

C 2=1000 
ASNM IN- 15000. 

C— 1 

INPUT  FORMAT 

15  WRITE  (5.  16) 

16  FORMAT  <///'  CODES  FOR  SELECTING  APPR.  PROB.  DIST.  '// 
1 1 5X»  'BINOMIAL',  12X.  '«1 

2/15X,  'POISSON',  13X.  '-2') 

READ  (5,*)  K 

IF  <  K. GT . 2. OR . K. LT. 1 )  GOTO  15 
22  WRITE (5.  21) 

21  FORMAT ( 10X,  'SELECT '/16X,  'SAMPLE  PLANS  ONLY  «1 ' 

1/16X, 'ASN  VALUES  ONLY  «2' 

2/16X, 'OR  BOTH  «3 ' ) 

READ  <5.  *>  KOPT 

IF  <  KOPT. QT. 3.  OR. KOPT.  LT. 1 )  GOTO  22 
WRITE  <5,  17) 

17  FORMAT < 10X,  'INPUT  ALPHA  ') 

READ  (5,*)  ALPHA 

WRITE  <5,  18) 

18  FORMAT < 10X,  'INPUT  BETA  ') 


JUC1)S7*MAIN 


0058 

0059 

0060 

0061 

0062 

0063 

0064 

0065 

0066 

0067 

0068 

0069 

0070 

0071 

0072 

0073 

0074 

0075 

0076 

0077 

0078 

0079 

0080 

0081 

0082 

0083 

0084 

0085 

0086 

0087 

0088 

0089 

0090 

0091 

0092 

0093 

0094 

0095 

0096 

0097 

0098 

0099 

0100 

0101 

0102 

0103 

0104 

0105 

0106 

0107 

0108 

0109 

0110 

0111 

0112 

0113 

0114 


READ  (5.*)  BETA 
WRITE  (  5<  51  ) 

51  FORMAT ( 10X. 'INPUT  PO  ') 

READ  <5.  *)  PO 

WRITE ( 5*  52) 

52  FORMAT ( 10X.  '  INPUT  PI  ') 

READ  (5.*)  PI 
WRITE!5. 58) 

58  FORMAT (  5X,  'INPUT  A  SEED  FOR  SINGLE  SAMPLING  NO.  '// 

1'  IF  NO  SEED  AVAILABLE  ENTER  ZERO  AS  THE  SEED  VALUE') 

READ! 5,*)  NS 
WRITE! 5# 59) 

59  FORMAT!  5X.  'INPUT  A  VALUE  FOR  !R1-C2>  '// 

1'  IF  R1-C2  THEN  THE  VALUE  WOULD  BE  0'// 

2'  IF  R 1 >C2  THEN  THE  VALUE  WOULD  BE  A  POSITIVE  NO.  '// 

3'  IF  R1!C2  THEN  THE  VALUE  WOULD  BE  A  NEGATIVE  NO. ') 

READ!  5.  *)  RU 

MC  1*=10.  0/ ! PI  /PO ) 

WRITE  ! 1. 53) 

53  FORMAT! ///10X.  'DEPT.  OF  ISE  ' 

1/.  10X. 'UNIVERSITY  OF  FLORIDA  '/ 

2/5X. 5! '*').  'DOUBLE  SAMPLING  SYSTEM '.  5! '*').  2X.  /  > 

WRITE  ! 5#  54 )  ALPHA.  BETA.  PO.  PI 
WRITE  !  1 .  54 )  ALPHA.  BETA.  PO.  PI 

54  FORMAT ! //10X*  'ALPHA  -'.F6.  4.  5X.  'BETA  »'»F6.  4. 

1/10X.  'PO  -'.F6.4.8X,  'PI  *'.F6.4) 

WRITE!5.  55)  RU 
WRITE!  1.  55)  RU 

55  FORMAT ! /5X»  'REJECTION  NO.  OF  FIRST  SAMPLE  !R1)  «  C2+! '.  13.  ')') 

5  C-C+l 

SINGLE  SAMPLING  PROCEDURE  BEGINS 

KK1*C+1 
10  NS-NS+1 

COMPUTATION  OF  LOWER  BOUND  OF  SINGLE  SAMPLING  PLAN 
IFiK.EQ.  1)  CALL  PR0BS1  INS.  PI.  C.  BXLEC.  N) 

IF!K. EQ.  2)  CALL  PR0BS2INS. PI.  C.  BXLEC.  N)  j 

IF ! BXLEC. CT.  BETA)  GOTO  10 
NLT-NS-5 
ML-MAXO! 1,  NLT) 

COMPUTATION  OF  UPPER  BOUND  OF  SINGLE  SAMPLING  PLAN 
20  NL-NL+1 

IF!K. EQ.  1 )  CALL  PR0BS1 <NL.  PO.  C, BXLEC) 

IF! K. EQ. 2)  CALL  PR0BS2!NL. PO. C. BXLEC ) 

IF!BXLEC.GE.  ! 1-ALPHA) )  GOTO  20 
NL-NL-1 

TEST  FOR  FEASIBILITY  j| 

IF! NS. GT. NL)  GOTO  5 


HfCAni  <r:  f .  V.'-  «P»'  TO8WB  V.  wjtj  <ru'  ViV*  AT-<fV 

S' 


^_r«”.  V*.  V--  w.  r.Tjv.-i-.  ¥--'*■<  i-j- 


K 


JGCU*J7*MALN 


01  1  5 
01  16 
0117 

one 

01  19 
0120 
0121 
0122 

0123  C 

0124  C 

0125  C 

0126 
0127 

0128  C 

0129 

0130  C 

0131 

0132 

0133  C 

0134  C 

0135  C 

0136 
0137 
0138 

0139  C 

0140  C 

0141 
0142 

0143  161 

0144 

0145 

0146  160 

0147  C 

0148 

0149 

0150 

0151  175 

0152  1 

0153  170 

0154  1 

0155 

0156  C 

0157  C 

0158  C 

0159 
0160 
0161 

0162  C 

0163  C 

0164  C 

0165 

0166 

0167 

0168 

0169  185 

0170 

0171  C 


'  )  ) 


WRITE  (5,90) 

WR I TE  (  1 , 90  ) 

90  FORMAT ( // 10X,  'SINGLE  SAMPLING  PLAN  '/'+'»9X.21( 

WR  ITE(  5,  91  >  C.NS,  NL 

WRITE (1,91)  C,NS,  NL 

91  FORMAT( 10X, 'ACCEPTANCE  NO. (C)  =  ',I2 

1. /10X,  'LOWER  BOUND  ON  N  (NS)  =  '.I4 

2,  /10X. 'UPPER  BOUND  ON  N  (NL)  =  '.I4) 

COMPUTATION  OF  DOUBLE  SAMPLING  PLAN  BEGINS:  FOR  EACH  VALUE  OF  C2 

IF ( C. LT. C2 )  MC=C+MC1-1 
C2=C 

R 1=C2+R 1 1 

DO  100  Kl  =  l, C2 
C 1 =K 1-1 

CALL  SUBROUTINE  TO  COMBUTE  THE  FIRST  SAMPLE  NUMBER 

CALL  TRYKNTRY,  Cl.  PI,  NS,  BETA.  K) 

N1=NTRY 

IF(NTRY. GT. NS)  GOTO  600 

WR  I  TE  <  5,  161) 

WRITE!  1. 161) 

FORMAT < /10X,  'DOUBLE  SAMPLING  PLANS'./) 

WRITE(5.  160)  C1.C2 
WRITE<  1.160)  C1.C2 

FORMAT (/10X.  'FOR  C1='.I2.  2X.  'C2='.  12. //) 

NTEMP=N1 

IF<KOPT.EQ. 1 )  WRITE! 5, 170) 

IF (KOPT.  EQ.  3)  WRITE<5.  175) 

FORMAT ( 10X,  '  <N1  )  '.  3X.  '  (N2S)  '.  3X.  '  ( N2L )  '.  4X, 

/) 

FORMAT <  1 IX,  '  (N1  )  '«  10X.  '  (N2S)  <  N2  < 

10X.  'PL'//) 

NTEMP1-NS 


'  ASNMAX  ',  5X, 
(N2L)  8X< 


'ASN ' 
'PS', 


COMPUTATION  OF  SECOND  SAMPLE  FOR  EACH  VALUE  OF  FIRST  SAMPLE 

ASN=FLOAT  <  NS ) * 1 0 
DO  190  IZ=NTEMP. NTEMP1 
I«IZ 

CALL  SUBROUTINE  TO  COMPUTE  SECOND  SAMPLE 

CALL  TRY2(NS.  NL.  K,  I.R1) 

IF (KOPT. NE.  1 )  GOTO  500 
WRITE  (5.185)  I,  N2S,  N2L,  PS,  PL 
WRITE!  1.185)  I,  N2S,  N2L,  PS,  PL 

FORMAT!  10X,  14.  13X,  14,  '  !  N2  <  ',  14,  4X.  F8.  6.  BX.  F8.  6) 

GOTO  190 
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g: 

[V 

LV 


i 


JQCDS7*MAIN 


•  0172 

•  0173 
0174 

.  0175 
0176 
0177 
0178 
0179 
0180 
0181 
0182 
0183 
0184 
0185 
0186 
0187 
0188 
0189 
0190 
0191 
0192 
0193 
0194 
0195 
0196 
0197 
0198 
0199 
0200 
0201 
0202 


C 

C 


C 

c 

c 


c 

c 


TEST  FOR  FEASIBILITY 

500  IF<N2S. GT.N2L)  GOTO  190 

IF < N2S. LT.  < C2-C 1 ) . OR.  I . LE. C2 ) 
ASNTEM=ASN 


GOTO  190 


C 

C 


CALL  SUBROUTINE  TO  COMPUTE  ASN(PO)  AND  ASNMAX  VALUES 

CALL  ASNN!MC,  NS. K.  I, KOPT, C1MIN. C2MIN.  N1MIN,  N2MIN.  ASNMIN) 
IFIKOPT.NE. 3)  GOTO  190 
WRITE!  5.  410)  I.  N2S.  N2L,  ASNMAX.  ASN 
WRITE!  1.  410)  I.  N2S.  N2L.  ASNMAX.  ASN 
410  FORMAT!  10X.  13.  4X»  13.  5X.  13.  5X.  2!F8.  4.  3X)  ) 

IF! ASN. GT.  ASNTEM)  GOTO  100 
190  CONTINUE 

100  CONTINUE 

600  IF ! C. LT. MC )  GOTO  5 

WRITE! 1.  601 )  ASNMIN, N1MIN, N2MIN. C1MIN, C2MIN 
WRITE! 5,  601)  ASNMIN,  N1MIN, N2MIN. C1M1N, C2MIN 

601  FORMAT!  10X,  'GLOBAL  MINIMUM  ASN!PO)« '. F8. 2. // 


110X, 'CORRESPONDING  N1 
21 OX,  'CORRESPONDING  N2S 
310X,  'CORRESPONDING  Cl 
410X. 'CORRESPONDING  C2 


STOP 

END 


15// 

15// 

12// 

12) 


'.\*.:-.-i 


4 

3 

5 

:j 

| 

\N 

ft 

li 


r.v 


is 


V 


w 

5 


0001 
0002 
0003 
0004 
0005 
0006 
0007 
0008 
0009 
0010 
001 1 
0012 
0013 
0014 
0015 
0016 
0017 
0018 
0019 
0020 
0021 
0022 
0023 
0024 
0025 
0026 
0027 


C 

C 

C 

c 

c 

c 

c 

c 

c 

c 


c 

c 


SUBROUTINE  TRY1 (NTRY,  Cl. P. NL. BETA,  K> 

THIS  SUBROUTINE  COMPUTES  FIRST  SAMPLE  NUMBER  OF  DOUBLE 
SAMPLING  PLAN  BY  AN  INTEGER  FORM  OF  BISECTION  METHOD 

INTEGER  Cl 


NLARGE=NL 

NSMALL=0 


NTR Y=  <  NSMALL+NLARGE ) / 2 . 0 

CALL  APPROPRIATE  PROBAILITY  SUBROUTINE  FOR  PROB. 
IFIK.EQ. 1)  CALL  PR0BS1 (NTRY, P, C 1 ,  BXLEC ) 

IF ( K. EQ. 2 )  CALL  PR0BS2(NTRY, P,  C 1 , BXLEC ) 

IF ( BXLEC. LE. BETA)  GOTO  50 
NSMALL=NTRY 
GOTO  25 
50  NLARGE=NTRY 

IF(NSMALL. NE. (NLARGE-1 > )  GOTO  5 
NTRY=NLARGE 
RETURN 
END 


10 


25 


CALCULATIONS 


M 


$4 


1 


9 


9 


0001 

0002 

0003 

0004 

0005 

0006 

0007 

0008 

0009 

0010 

0011 

0012 

0013 

0014 

0015 

0016 

0017 

0018 

0019 

0020 

0021 

0022 

0023 

0024 

0025 

0026 

0027 

0028 

0029 

0030 

0031 

0032 

0033 


SUBROUTINE  PR0BD1 (Nl, N2, P, DPROB. K. R 1 ) 

THIS  SUBROUTINE  COMPUTES  DOUBLE  PROBABILITIES  FOR 
COMPUTING  SECOND  SAMPLE  NUMBER  OF  DOUBLE  SAMPLING 

C0MM0N/BLK6/C 1 » C2 
INTEGER  C1.C2,  R1 


NUMBER 


IF ( K. EG.  1 )  CALL  PR0BS1 < Nl , P, C 1 , BXLEC ) 

IF (K. EG. 2)  CALL  PR0BS2(N1,  P, Cl. BXLEC ) 

DPROB=BXLEC 

TEMP=BXLEC 

NTEMP=C 1+1 

KTEMP=R1-1 

DO  10  IX=NTEMP»  KTEMP 
I  =  IX 
J=C2— I 

IF(K. EG. 1)  CALL  PR0BS1 (Nl. P, I.  BXLEC) 
IF<K. EQ. 2)  CALL  PR0BS2(N1» P, I, BXLEC) 
PROB 1=BXLEC-TEMP 
TEMP=BXLEC 

IF(K. EG. 1 )  CALL  PR0BS1 ( N2i P, J, BXLEC ) 
IFIK.EG.  2)  CALL  PR0BS2 ( N2.  P,  J,  BXLEC  ) 
DPR0B«DPR0B+(PR0B1*BXLEC) 

10  CONTINUE 

RETURN 

END 


SUBROUTINE  TRY2<NS»  NL*  K,  J#  R1 ) 


THIS  SUBROUTINE  COMPUTES  THE  SECOND  SAMPLE  NUMBER  OF 
THE  DOUBLE  SAMPLING  NUMBER  BY  AN  INTEGER  BISECTION 
METHOD.  SEVERAL  TESTS  ARE  DONE  TO  LOCATE  THE  PARAMETER 
AT  ITS  TRUE  POSITION. 

INTEGER  Cl, C2» R1 

COMMON/BLK 1 /N2Sj N2L 
C0MM0N/BLK2/PS.  PL 
C0MM0N/BLK3/N1 
C0MM0N/BLK4/ALPHA.  BETA 
COMMON/BLK5/PO,  PI 
C0MM0N/BLK6/C1,  C2 

K 1 =C 1+1 

SET  LIMITS  FOR  COMPUTING  N2S 

NSMALL=NS-J 

NLARGE=NSMALL 

INDEXING  TO  SPECIFY  WHAT  BOUND  (N2S  OR  N2L)  IS  BEING 
COMPUTED 


INITIAL  TEST  AT  EACH  LIMIT 

CALL  PROBD1 <  J,  NSMALLj PI,  DPROB, K. Rl) 
IF(DPROB.LE. BETA)  GOTO  55 

BISECTION  METHOD 

NLARGE=NL 

5  NTRY=*<  NSM ALL+NLARGE )  /2 .  O 
GOTO  <10,  20).  I 

10  CALL  PROBDKJ,  NTRY,  PI,  DPROB,  K,  Rl) 

IF < DPROB. LE. BETA)  GOTO  50 
GOTO  15 

20  CALL  PR0BD1<J, NTRY, PO, DPROB, K. Rl ) 
IF(DPROB.LT.  <1 -ALPHA ) )  GOTO  50 
15  NSMALL*=NTR  Y 
GOTO  25 

50  NLARGE=NTRY 

25  I F  <  <  NLARGE-NSMALL ) . GT. 1 )  GOTO  5 

CHECK  THE  INDEX  TO  FIND  WHERE  THE  PROCESS  IS 
GOTO  <55,  60),  I 


7 


.VftVi 


V»V: 


FRY2 


0058 

0059 

0060 

0061 

0062 

0063 

0064 

0065 

0066 

0067 

0068 

0069 

0070 

0071 

0072 

0073 

0074 

0075 

0076 

0077 

0078 

0079 

0080 

0081 


CHANGE  THE  INDEX  AFTER  N2S  COMPUTATION 
55  1*1+1 

TESTING  EACH  POSSIBLE  CASES  TO  LOCATE 
THE  LOWER  BOUND  AT  ITS  TRUE  POSITION 

N2S=MAX0<0» NLARGE ) 

CALL  PR0BD1 ( J, NLARGE. PI. DPROB.  K.  R 1 ) 

PS*DPROB 

MTEMP=NLARGE-5 

NSMALL=MAXQ ( 0 .  MTEMP ) 

NLARGE=NL 
GOTO  5 

60  N2L=N5MALL 

CALL  PR0BD1  ( J,  NSMALL.  PO.  DPROB,  K,  R 1 ) 
PL=DPROB 

CALL  PR0BD1 ( J. NLARGE. PO. DPROB. K. R1 ) 

IF ( DPROB. GE. < 1 -ALPHA ) )  N2L=NLARGE 
IF< DPROB. GE.  < 1 -ALPHA > >  PL*DPROB 


110  RETURN 
END 


,r*,i 


v.v 


0001 

0002 

0003 

0004 

0005 

0006 

0007 

0008 

0009 

0010 

OOll 

0012 

0013 

0014 

0015 

0016 

0017 

0018 

0019 

0020 

0021 

0022 

0023 

0024 

0025 

0026 

0027 

0028 

0029 

0030 

0031 

0032 

0033 

0034 

0035 

0036 

0037 

0038 

0039 

0040 

0041 

0042 

0043 

0044 

0045 

0046 

0047 

0048 

0049 

0050 

0051 

0052 

0053 

0054 

0055 

0056 

0057 


C 

c 

c 

c 

SUBROUTINE  PR0BS1 (NN« P, C» BXLEC ) 

C 

C  THIS  SUBROUTINE  COMPUTES  CUMULATIVE  BINOMIAL 

C  PROBABILITIES 

C 

INTEGER  C 

DOUBLE  PRECISION  SUMLOG 
CQMM0N/BLK7/SUMLQG ( 1500) 

C0MM0N/BLK8/N 

C 

C 

Q=1 . -P 
C 

C  BINOMIAL  PROB.  WHEN  C=0 

C 

CSUMS=Q**NN 

C  WRITE (6,  500)  CSUMS 

IF  <C.EQ. 0)  GOTO  333 
C 

C  AVOID  RECOMPUTING  SUMLOG ( I ) 'S  ALREADY  IN  MEMORY 

C 

IF  (N-NN)  100i 21 1«  211 
100  M=N+1 
C 

C  COMPUTE  N  SUMLOGS-EGU I VALENT  TO  N— FACTORIAL 

C 

IF  <M.GT.  1)  GOTO  110 
SUMLOG ( 1 )=0. 

IF(NN. LE. 1 )  GOTO  211 
M=2 

110  DO  111  I=M»  NN 

SUMLOG ( I ) =DLOG 1 O ( DFLOAT ( I ) ) +SUMLOG ( I - 1 ) 

111  CONTINUE 
C 

C  COMPUTE  C  SUMS-EQU I VALENT  TO  SSUM  OF  PROB.  COMP IN. 

C  I.E.  CUMULATIVE  BINOMIAL  DISTRIBUTION  COMPUTATION 

211  IF(NN. GT. N)  N=NN 
C 

C  DETERMINE  BEST  NUMBER  HANDLING  LOOP 

C 

IF  (NN. GT. 300)  GOTO  300 
DO  222  K=1 »  C 

CSUMS=10.  **< SUMLOG < NN > -SUMLOG ( NN-K > -SUMLOG < K > ) 

1  *P**K*Q**<NN-K)+CSUMS 

222  CONTINUE 

C  WRITE<6,  501)  CSUMS 

C  501  FORMAT < 5X«  'XXX  '.FB.  6) 

GOTO  333 
C 

C  LOOP  FOR  LARGE  EXPONENTS 

C 

300  DO  322  K=l»  C 
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PR0BS1 

00S8 

C  SUMS= 1 0. ** ( SUMLOG ( NN ) -SUMLOG ( NN-K ) -SUMLOG  <  K ) 

0059 

1  +K*DL0G10<DBLE(P) )  +  <  NN-K  >  *DL0G 1 0  <  DBLE ( G ) ) )+CSUMS 

0060 

C 

WRITE (6, 501)  CSUMS 

0061 

C 

500 

FORMAT ( 10X»  F8.  6 ) 

0062 

322 

CONTINUE 

0063 

0064 

C 

333 

BXLEC  =  CSUMS 

0065 

RETURN 

0066 

END 

v 
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0001 

C 

0002 

c 

0003 

c 

0004 

c 

0005 

0006 

r 

SUBROUTINE  PRQBS2<NN»  P.  C»  BXLEC  > 

0007 

c 

THIS  SUBROUTINE  COMPUTES  CUMULATIVE  POISON 

0008 

c 

PROBABILITIES 

0009 

c 

0010 

INTEGER  C 

OOll 

pp=»p#nn 

0012 

TERM=1 . 0 

0013 

SUM=TERM 

0014 

c 

0015 

IFIC.EQ. 0)  GOTO  110 

0016 

DO  100  1=1. C 

0017 

TERM=TERM*PP/I 

0018 

SUM=SUM+TERM 

0019 

100 

CONTINUE 

0020 

c 

0021 

1  io 

BXLEC=SUM/EXP(PP ) 

0022 

c 

0023 

RETURN 

0024 

END 
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SUBROUTINE  ASNN ( MC. ,  NS,  K,  N1  1  ,  KOPT  ,  C1MIN,  C2M IN.  N1MIN,  N2MIN.  ASNMIN) 

THIS  SUBROUTINE  COMPUTES  ASN(PO)  VALUES  AND  ASNMAX 
VALUES. 

DOUBLE  PRECISION  SUMLOG 
INTEGER  C1MIN, C2MIN 
C0MM0N/BLK1/N2S,  N2L 
COMMON/BLK3/N1 
C0MM0N/BLK4/ALPHA,  BETA 
C0MM0N/BLK5/P0,  PI 
C0MM0N/BLK6/U.  12 
C0MM0N/BLK7/SUML0G ( 2500 ) 

COMMON/BLKB/N 
C0MM0N/BLK11/ASN,  ASNMAX 

INITIALIZATION 

COMPUTE  P*  (MAXIMUM  PROB.  FOR  ASNMAX) 

J=I1+1 
XXX=0. O 

IF(Il.GT.O)  X  X  X=SUMLOG (II) 

AKONST=10. **( SUMLOG ( 12 ) +SUMLOG ( N1 1—12—1 )-XXX-SUMLOG(Nl 1-1 1-1 ) ) 
TEMP= 1 . O/FLOAT ( 1 2- 1 1 ) 

AKONST =AKONST  **TEMP 
PSTAR=AKONST / ( 1.  O+AKONST) 

IF ( K. EQ. 1 )  CALL  PR0BS1 (Nil/ PSTAR.  12.  BXLEC ) 

IF(K.EQ.2)  CALL  PR0BS2(N1 1,  PSTAR,  12. BXLEC ) 

TEMP=BXLEC 

IF (K. EQ. 1 )  CALL  PR0BS1 (N1 1 , PSTAR. I 1, BXLEC ) 

IF(K. EG. 2)  CALL  PR0BS2 (N1 1 . PSTAR.  1 1 , BXLEC ) 

TEMP 1=TEMP-BXLEC 
ASNMAX=FL0AT(N11 )+N2S*TEMPl 

IF(K. EQ. 1)  CALL  PR0BS1 (N1 1 , PO. 12. BXLEC ) 

IF (K. EQ. 2)  CALL  PR0BS2(N1 1,  PO,  12, BXLEC ) 

TEMP=BXLEC 

IF(K. EQ. 1 )  CALL  PR0BS1 (N1 1 , PO. I 1 . BXLEC ) 

IF(K. EQ. 2)  CALL  PR0BS2(N1 1.  PO.  1 1. BXLEC ) 

TEMP 2=TEMP -BXLEC 
ASN=FL0AT(N1 1 )+N2S*TEMP2 

IF( ASNMAX. GT.  NS.  OR.  ASN.  GT.  ASNMIN)  GOTO  100 

ASNM I N= ASN 

N1MIN=N1 1 

N2MIN-N2S 

C1MIN*I1 

C2MIN*I2 

100  IF ( KOPT. NE. 2 )  GOTO  500 

WRITE (5.  110)  Nil.  N2S,  TEMP  1 ,  ASNMAX.  TEMP2.  ASN 
110  FORMAT  ( /  /10X,  2(13.  3X),  2(F6.  4.  3X,  F8.  4,  3X)  ) 


8  * 


C 

C 


